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Abstract 

Low frequency noise in current biased Lao.82Cao.isMn03 single crystals has been investigated in 
a wide temperature range from 79 K to 290 K. Despite pronounced changes in magnetic properties 
and dissipation mechanisms of the sample with changing temperature, the noise spectra were found 
to be always of the 1/f type and their intensity (except the lowest temperature studied) scaled as a 
square of the bias. At liquid nitrogen temperatures and under bias exceeding some threshold value, 
the behavior of the noise deviates from the quasi-equilibrium modulation noise and starts to depend 
in a non monotonic way on bias. It has been verified that the observed noise obeys Dutta and Horn 
model of 1/f noise in solids. The appearance of nonequilibrium 1/f noise and its dependence on 
bias have been associated with changes in the distribution of activation energies in the underlying 
energy landscape. These changes have been correlated with bias induced changes in the intrinsic 
tunneling mechanism dominating dissipation in Lao.82Cao.isMn03 at low temperatures. 

PACS numbers: 72.70.+m Noise processes and phenomena, 72.15.-v Electronic conduction in metals and 
alloys, 73.43.Jn Tunneling, 75.47.Gk Colossal magnetoresistance 
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I. INTRODUCTION 



Ultraslow relaxation in a variety of physical and nonphysical systems, manifesting itself 
in 1/ f a noise power spectral density (PSD) of the fluctuating observable quantity, with 
a ~ 1, continues to be one of the most exciting and puzzling effects attracting the attention 
of many research groups.-^ The most commonly investigated low frequency fluctuations in 
solid state system are conductivity fluctuations. It is generally recognized that 1/ / conduc- 
tivity noise in solids arises from superposition of many elementary Lorentzian contributions 
of individual fluctuators with well defined and specifically distributed characteristics relax- 
ation rates.- Theoretical models of 1/f noise in solids relate changes in the noise PSD to 
changes in the energy landscape in which charge carriers are moving and dissipating energy. 
Therefore, studies of noise in solid state systems are important not only from the point of 
view of improving the performance of solid state devices but also provide a unique tool for 
understanding the nature and dynamics of transport processes in an investigated system.- 

The issue of noise in colossal magnetoresistance (CMR) manganites with general formulae 
Ri-^A^MnOs, where R is a rare-earth ion and A a divalent or tetravalent cation, was ad- 
dressed in many publications in the last decade. Almost all of them revealed prominent broad 
band conductivity fluctuations with PSD following l/f a law. Low-doped Lai-^Ca^MnOs 
(LCMO) system exhibits metal-to-insulator (M-I) transition in the vicinity of Curie tem- 
perature Tq. Above the M-I transition the conductivity of LCMO is dominated by hop- 
ping mechanism which strongly depends on oxygen stoichiometry and the tolerance factor. 
Sources of 1/f noise at Tq were therefore initially associated with dynamics of oxygen va- 
cancies.- 

Early experiments concentrated on the noise peak associated with M-I transition, in the 
vicinity of which a 3 to 4 order of magnitude increase of the noise level has been observed 
in Ca and Sr doped manganites, see e.g.- - -. The noise peak around M-I transition has 
been interpreted in terms of percolative nature of the transition between charge ordered 
insulating and metallic ferromagnetic states, 6 or in terms of magnetic fluctuations coupled 
to the resistivity- It became however soon evident that this unusually large 1/f noise is not 
an intrinsic property of doped manganites. It was found that 1/f noise peak appears only 
in notably strained films while it is completely absent in single crystals and in high quality, 
almost strain-free, epitaxial films.- - — Small noise peak showing sometimes around Tq in high 
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quality samples could be easily suppressed by application of a small magnetic field and was 
attributed to magnetic fluctuations in the vicinity of the phase transition. Experimental 
observations of direct correlation between magnetic noise and magnetoresistance allowed to 
conclude that the enhanced 1/f noise originates in charge carriers density fluctuations.— 
Current, temperature, and magnetic field dependence of noise in grain boundary junctions, 
provided clear evidence that 1/f noise is caused by localized states with fluctuating magnetic 
moments in heavily disordered grain boundary regions.— 

Conductivity noise with 1/f spectrum is generally related to resistance fluctuations which 
are measured by applying dc current and recorded as voltage fluctuations. When the resis- 
tance fluctuations are just probed by current, and not influenced by its flow, then PSD of the 
noise scales as the square of the bias current. Such modulation noise is referred to as quasi- 
equilibrium 1/f noise. There is a mounting experimental evidence that quasi-equilibrium 
1/f noise in doped manganites is accompanied by the noise that arises from, or is directly 
modified by, the passage of current through the sample.- 1 ^ 1 ^ This type of noise is referred 
to as nonequilibrium 1/f noise, and its dependence on bias is considerably different from a 
quadratic one. 

Our system of interest is Lao.82Ca .i8Mn0 3 (0.18 LCMO) single crystal, a highly corre- 
lated electron system with a complex interplay of charge, orbital, and spin ordering, what 
leads to remarkable changes in transport characteristics and magnetic state of the sample 
with decreasing temperature. In this way 0.18 LCMO compound provides us with a unique 
opportunity of studying noise properties in the same solid-state sample, in markedly dif- 
ferent transport regimes and magnetic properties, using temperature and current bias as 
independent factors modifying the state of the active fluctuators. 

In this paper we demonstrate that equilibrium 1/f noise in Lao^Cao.isMnOs single 
crystals preserves its character under changing temperature, despite significant changes in 
magnetic and transport properties of the system. Bias dependent nonequilibrium 1/f fluc- 
tuations appear only at low temperatures, well below the Curie temperature. Properties of 
the nonequilibrium 1/f noise differ significantly from those of quasi-equilibrium one. A com- 
bined analysis of noise and transport characteristics permits us to associate marked changes 
in the noise behavior to changes in the low temperature intrinsic tunneling mechanism. 
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II. EXPERIMENTAL BACKGROUND 



Lao^Cao.isMnOs crystals were grown by a floating zone method using radiative heating. 16 
The crystallographic orientations of the crystal and wafers were determined by Laue method 
with accuracy of ±2°. X-ray data of the crystal were compatible with the perovskite struc- 
ture orthorhombic unit cell, a = 5.5062 A, b = 7.7774 A, c = 5.514 A. The as-grown crystal, 
in form of a cylinder, about 4 cm long and 4 mm in diameter, was cut into individual small 
rectangular 6x3x2 mm 3 bars, with the longest dimension along the < 110 > crys- 
tallographic direction. Current and voltage leads were indium soldered to gold/chromium 
contacts deposited by thermal evaporation in vacuum. 

For noise measurements the sample was thermally anchored to the sample holder of 
a variable temperature liquid nitrogen cryostat. Conductivity noise was measured in a 
conventional 4-point contact arrangement by biasing the sample with dc current supplied by 
high output impedance current source and measuring the resulting voltage fluctuations. Four 
in-line contacts were placed along the longest dimension of the bar. The separation between 
voltage contacts was 0.3 mm. The voltage signal was amplified by a home made room 
temperature low noise preamplifier, located at the top of the cryostat, and further processed 
by a computer assisted digital signal analyzer. To eliminate environmental interferences and 
noise contributed by the measuring chain, the PSD measured at zero current was subtracted 
from the data obtained at a given current flow for each measurement. 

Remarkable changes in transport properties and magnetic state of 0.18 LCMO system 
upon changing temperature are illustrated in Fig. [TJ At high temperatures the sample 
is in paramagnetic insulating state and the resistivity, dominated by hopping mechanism, 
increases with decreasing temperature. The resistivity reaches a pronounced maximum re- 
lated to the metal-insulator (M-I) transition at T = Tui = 180 K. The maximum appears at 
temperature very close to the Curie temperature Tq of the paramagnetic-to-ferromagnetic 
transition. Tq determined by independent magnetization measurements is 180 ±1 K. Intrin- 
sic phase separation (PS) associated with metal-insulator (M-I) transition at T ~ Tq leads to 
percolation conductivity in the ferromagnetic (FM) state at T < Tq. Therefore, at temper- 
atures below the magnetic ordering temperature the resistivity decreases in a metallic-like 
way with dR/dT > 0, although its absolute value is much higher than that of common met- 
als. With further temperature decrease, the resistivity reaches a shallow minimum around 
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Temperature [K] 

FIG. 1: Color online. Temperature dependence of the resistance of Lao.82Cao.isMn03 crystal at 
zero applied magnetic field. The room temperature resistivity of the crystal is 2.4 Qcm. 

T ~ 120 K, followed by a strong upturn at temperatures below some 100 K. 

Detailed magnetic measurements of our sample did not reveal any peculiarities at tem- 
peratures below Tc, suggesting that in the investigated temperature range there is only one 
magnetic transition at T = Tc- Low temperature resistivity is dominated by tunneling 
across intrinsic barriers associated with extended structural defects, such as twins and grain 
boundaries, and/or with inclusions of insulating FM phase interrupting metallic percolating 
paths.— 

Despite changes in the dissipation mechanism and magnetic state of the sample, the 
experimentally observed voltage noise was found to have 1/f spectrum in the entire inves- 
tigated temperature and bias range. To get a first insight into the nature of the observed 
noise, we have measured the dependence of the noise intensity on bias current. We found 
that for all temperatures, with the exception of the 79 K data, for currents exceeding / ~1 
mA, the noise intensity scales as the square of the bias current, as illustrated in Fig. |2j 
Proportionality of 1/f voltage noise to the current squared indicates that the noise is due 
to bias independent resistivity fluctuations which are unaffected by the current flow. 
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FIG. 2: Color online. Current dependence of the noise intensity at / = 1 Hz as recorded at various 
temperatures. 

At low temperatures, where electrical transport in our sample is dominated by tunneling 
mechanism, 17 at currents exceeding some threshold current / ~ 1 mA, the proportionality of 
the voltage noise to I 2 breaks down and the resistivity fluctuations start to be influenced by 
the current. This is a manifestation of a nonequilibrium 1/ / noise similar to that frequently 
observed in nonlinear, nonohmic systems.- 8 However, in a difference to other nonlinear 
systems, the crossover to nonequilibrium noise in our sample is not associated with the 
onset of strong nonlinearity in the current-voltage (I — V) characteristics.— Moreover, the 
intensity of nonequilibrium noise initially decreases with increasing bias, to change in a 
nonmonotonic way with further increase of the bias. 

The puzzling properties of the low temperature noise in Lao^Cao.isMnOs single crys- 
talline sample should be related to the underlying changes in its transport properties. An 
efficient way to approach the problem is to detect possible changes in the energy environ- 
ment of the fluctuators responsible for the noise. A tool for such analysis is provided by 
relevant theoretical models of nonexponential relaxation kinetics resulting in 1/ f noise. 1,2 
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III. THEORETICAL BACKGROUND 



A. 1/f noise in solids 

1/f noise spectra in condensed-matter systems are believed to come from assembly of 
elementary fluctuators with well defined characteristic relaxation rates A. The spectrum of 
each elementary fluctuator is Lorentzian, and the resulting normalized PSD of the fluctuating 
quantity X, 



will have l/f a form, with a = 1 for the distribution function P(X) oc 1/A. If the elementary 
fluctuators are thermally activated, the characteristic rate A(A) = Aoexp(— A/k B T), and 
the required distribution of A is provided by flat distribution of activation energies, -P(A) = 
const. 

Dutta and Horn, (DH) have shown that l/f a spectra with a ~ 1 arise not only for 
P(A) = const, but also for the distribution function P(A) that does not vary much in the 
range of A from |A;BTln(Ao//i)| till \k B T\n(X / f 2 )\, where fi and f% are the lowest and the 
highest frequency of the measurements 1 . For a detailed discussion of DH model assumptions 
see reference Q 

When P(A) is a slowly varying function of the activation energy A then 



In the majority of solid state systems Ao is a constant of the order of phonon frequencies 
Ao ~ 10 12 — 10 13 Hz.- But if a disordered system undergoes a phase transition at transition 
temperature Tq, then for the noise spectra in the vicinity of Tq we obtain the same Eq. [2] 
but with k B T replaced by l/£, e , where £ is the correlation length and 9 is a dynamic critical 
exponent.-^ As a result, while analyzing experimental data using Eq. [2] one may obtain 
unreasonably large values of Ao- 

In the experiments we have measured the noise in the frequency range from /o = 1 
Hz to / = 10 2 Hz. Therefore, ln(///o) < 5 may be considered as a small addition to 
ln(Ao//o) ~ 25 — 30 and we can further expand the P(A = k B T\n(\o/ f)) function from the 
right-hand side of Eq. (J2J) as 




(1) 



S(f,T)f oc k B TP(A = k B T ln(A //)). 



(2) 



P(k B T\n(X /f)) « P(fc B Tln(A //o)) 



dP 



k B Tln(f/f ). 



(3) 



dA A=k B Tln{X /fo) 
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Within the DH approximation, at low frequencies, S(f,T)f-vs.-\ia(f/f ) dependence should 
be linear at a fixed temperature. It follows from Eqs. [2]and[3]that the slope of S(f, T)f-vs- 
ln(/// ) plot is proportional to dP(A)/dA. This slope is directly related to PSD exponent 
a through a general relation 

,, T1 , 1 d(S(f,T)f) 

DH model leads to a specific relation between temperature derivative of the spectral 
density dln(S(f,T)f)/dlnT and the derivative d\n(S(f,T)f)/d\nf 

dln(S(f,T)f) /f , d\R(S(f,T)f) 

— 1 - ln(A //) — . (5) 

Equation (JSJ) contains crucial reciprocity between the frequency and temperature depen- 
dence of the noise magnitude and is frequently used as a self-consistency test for the validity 
of DH approach for a given physical system. The relevant test consists in verifying whether 
the dependence of the noise intensity on frequency at a constant temperature, and on tem- 
perature at a constant frequency, are consistent with each other. If this is the case, one 
may use Eq. (j2J) to attribute a single distribution of activation energies to the ensemble of 
active fluctuators. Flat distribution of temperature independent activation energies in DH 
model gives rise to a pure 1/ / spectrum and a linear temperature dependence of the noise 
level at low frequencies. Departures of a from unity indicate non-zero derivatives of energy 
distribution dP(A)/dA, and excess of high (a < 1), or low energy (a > 1) fluctuators in 
the ensemble. 



B. Inelastic transport across an inhomogeneous medium. 

At lowest temperatures studied in the experiment, charges in the investigated system 
are transmitted by tunneling across classically forbidden regions. Electron wave functions 
in tunneling barriers decay exponentially over the characteristic length I. When the decay 
length / exceeds the barrier thickness d, then direct tunneling, both elastic and inelastic, 
dominates the conduction mechanism. Elastic contribution to the differential conductance 
G(V) = dI(V)/dV- versus- voltage characteristic is a parabola 20 , with 

G cl (V) = G + const • (V - V ) 2 , (6) 

where G is proportional to exp(— 2d/t). 
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The inelastic contribution to the differential conductance is given by Gi n (V) oc 

eV 

J a 2 F(u)du with the Eliashberg electron-boson interaction function a 2 F(u).— For most 
o 

metals, in the case of phonons as mediating bosons, the shape of the average Eliashberg 
spectral function resembles that of the phonon density of states F(u).— Let us look, for 
example, at F{uS) obtained in neutron scattering experiments for a polycrystalline sample 
of Lao.625Cao.375Mn03.— One can see that, within the range of 10 to 70 meV, the average 
deviation of the phonon density of states from its mean value does not exceed 20 %. There- 
fore, in this energy interval one may use the approximation F(u) = const, which results in 
a linear behavior of G in vs. V; G in (V) oc \ V\.— Such situation has been frequently observed 
in disordered perovskite oxides, high Tq cuprates and various CMR manganites. — 

With increasing thickness of the tunneling barrier, hopping along chains of localized states 
is favored. While hopping along localized states path, a carrier does not cross quantum- 
mechanically the entire distance between the electrodes, but rather jumps from the junction 
electrode to the first state, lose the phase memory, moves to the second nano-island and, 
eventually, after completing all the hopping path, jumps to the opposite electrode. Such 
case was considered by Glazman and Matveev^ 1 ^ in their model (GM) of indirect tunneling 
in disordered materials. The GM model applies very well to low temperature transport in 
perovskite manganites. iL2Zi2& 

In the GM model, temperature and voltage dependencies of the tunnel conductance are 
expressed as multistep tunneling via N localized states: 

oo 

G(V) = G +Y, G n(V,T), (7) 

N=l 

where conductance G represents bias and temperature independent elastic tunneling term, 
while Gn describe tunneling through N > 1 localized states. 

G N (V) = a N V {N -^ ] for eV > k B T, (8) 
G N (T) = b N T {N -^ for eV < k B T } (9) 

where coefficients a at and depend exponentially on barrier thickness. 

Resuming, the differential conductance of a thin, inhomogeneous, insulating barrier is a 



power function of the voltage, see also Ref. 
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G in (V) = G + const • V n , (10) 
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where index n characterizes the tunneling regime: n = 2 corresponds to elastic tunneling 
with the energy relaxation in the conducting regions of the system, whereas other n are 
signatures of inelastic tunneling in which an electron losses its energy inside the insulating 
region. By finding the value of index n from experimental data one can infer information 
about the physics of electron transport processes across the dielectric layer. 

IV. RESULTS AND DISCUSSION 
A. Applicability of DH model 

It should be noticed that deviations or inconsistencies with DH model in solid state 
systems are rare.- In particular, it has been already demonstrated that noise in various half- 
metallic oxides obeys the DH model.- Nevertheless, as a first step we perform tests of the 
applicability of DH model to our data by checking the reciprocity condition expressed in 
Eq. (jSJ) which relates the temperature and frequency derivatives of the dimensionless PSD 
at fixed T and / values. Figure |3] compares left and right sides of Eq. (jSJ) for frequency 
equal to 2 Hz, current bias of 1.5 mA and all temperatures studied. The overall agreement 
is good, including also the 79 K data, and even the change of the sign around 175 K is well 
reproduced in the both dependencies. 

Canonical DH approach assumes that P(A) is temperature independent. We find a 
small violation of this assumption for temperatures around Curie temperature Tc- The 
parameter Ao estimated from the best agreement, in the least-squares sense, between two 
dependencies plotted in Fig. [3] was found to be Ao = 4.3 x 10 14 Hz. Thus obtained Ao is 
physically reasonable, although it is about two orders of magnitude higher than a typical 
phonon frequency. We attribute this discrepancy to the existence of the phase transition 
in the investigated system at T = 180 K.-^ Indeed, separate estimations of Ao for each 
individual temperature render values close to the expected phonon frequency at high and 
low temperatures, while values of Ao for temperatures close to Tc are unphysical and diverge. 

We conclude that the frequency-temperature dependence reciprocity condition is fulfilled 
and that the observed noise does obey the DH model in entire investigated temperature 
range. We are allowed therefore, to assume that the observed noise is due to thermally 
activated kinetics of active fluctuators with the distribution of activation energies P(A) given 
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FIG. 3: Color online. Verification of the DH temperature-frequency reciprocity. Comparison of 
the left- and right-hand sides of Eq. ([5]) at / = 2 Hz and I = 1.5 mA. 

by Eq. (J2J), and the derivative of activation energy distribution dP(A)/dA proportional to 
the slope of S(f,T)f/V 2 -vs.-hi(f/f ) characteristics. 

B. Quasi-equilibrium noise. 

Figure H] shows the data from Fig. [2]replotted as a function of bias voltage. Within the 
experimental accuracy data recorded at all temperatures, with the exception of T = 79 K, 
collapse into single line Sy oc V 2 . This is equivalent to the statement that Sy/V 2 = const, 
meaning that the reduced spectral density of resistance fluctuation S r = Sr/R 2 = Sy/V 2 
is temperature independent. Temperature independent S r implies that the dependence of 
resistance fluctuations 5R on temperature follows the temperature dependence of the resis- 
tance of the sample R. Experimentally revealed condition 5R oc R assures as well that the 
coupling of the fluctuators to resistance fluctuations is temperature independent. Observe 
that temperature independence of the coupling of system fluctuations to the measured phys- 
ical parameter is one of the most problematic condition assumed in the Dutta-Horn model 
of 1/ f noise and its violation frequently leads to some deviations of the experimental data 
from the simplest version of the model.- On the other hand the temperature independence 
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FIG. 4: Color online. Voltage dependence of the noise intensity at / = 1 Hz as recorded at various 
temperatures. 

of the above coupling is quite surprising result if one takes into account pronounced changes 
of the magnetic state and dissipation mechanism in the sample with changing temperature. 

Upper panel of Figure El shows S(f, T)f data recorded at 150 K for different bias currents 
and plotted as a function of ln(///o), while the lower panel shows the same data normalized 
by dc voltage squared, S(f,T)f/V 2 . First of all, one observes that Sy(f,T)f vs. ln/// at 
a fixed temperature is indeed a linear function as predicted by Eq. ([3]) (similar results were 
obtained for other currents as well). Moreover, all normalized spectra recorded at different 
currents collapse to a single line, consistently with the I 2 dependence of the equilibrium 
noise. The conductivity noise at 150 K is thus only probed and not modified by current 
flow. By fitting the normalized data to a linear function (a solid straight line in Fig. [5]) one 
obtains the derivative d(S(f, T)f )/dln f. Therefore, we may conclude immediately that 
at T = 150 K bias current / has no influence whatsoever on the distribution of activation 
energies P(E) since the slope of each spectrum is the same within the experimental accuracy, 
meaning that dP(A)/dA does not depend on bias. The same independence of the slope 
factor, and consequently P (A) function, on current has been confirmed for all temperatures 
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FIG. 5: Color online. Upper panel: Example of a linear dependence of the product Sy(f,T)f 
as a function of ln///o recorded for various currents at T = 150 K. Lower panel: Data from 
the upper panel after normalization of the noise intensity by V 2 . Observe that all data from the 
upper panel collapse to a single line after the noise normalization. Solid straight line is a linear-fit 
approximation of the experimental data for all current biases. 
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FIG. 6: Color online. Effect of the current bias on the dimensionless noise intensity recorded at 79 
K. 

with exception of 79 K data. Here, even a brief look at the plot S(f, T)f /V 2 -vs-ln f from 
Fig. Improves that at 79 K, in a marked difference to higher temperatures, the distribution 
function and its derivative are strongly influenced by the current flow. This is reflected in a 
strong dependence of the spectral exponent a on current flow at 79 K, see Eq. HI in contrast 
to a very weak current influence on a at higher temperatures. 



C. Nonequilibrium noise. 

In order to measure resistance fluctuations in the equilibrium state, a small current is 
applied to a sample with the only aim to convert the resistivity fluctuations into measurable 
voltage noise. In this case the dimensionless PSD spectrum Sy(f,T)f/V 2 is a constant 
which does not depend on voltage, or current, bias. 

It can be clearly seen in Figs. [2] and 0J that 1/f noise at 79 K looses its equilibrium 
character at bias exceeding 1 mA, or equivalently, 200 mV. Above the threshold bias the 
noise is no more proportional to bias squared. The nonequilibrium noise intensity changes 
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in a nonmonotonic way with increasing bias and decreases with increasing bias at certain 
bias ranges. 

There are several reasons why it may happen. We shall discuss here four possibilities 
which can have direct relation to our study. The first two are artifacts associated with 
shunting parts of the sample under the voltage contacts by the contact pads or simple 
overheating of a sample due to dissipation at high current biases. The third possibility is 
a direct influence of the external bias on elementary fluctuators. Let us consider a single 
two-state system with equal activation energies of both states, symmetric two well energy 
structure, at zero bias. If the energy structure of the fluctuator is stressed by the bias, 
the energy difference between metastable states 1 and 2 will increase with increasing bias. 
The total switching rate is a sum of the two switching rates A ii2 = X exp(—A 12 /k B T) for 
two metastable states.— This mechanism may lead to increasing asymmetry in the energy 
structure of the involved fluctuators and, consequently, in a decrease of noise intensity 
with increasing bias. A fourth possibility of remarkable bias enforced modifications of the 
noise PSD, which would be most important for inhomogeneous samples, are changes in the 
mechanisms of electron transmission through the investigated sample. 

Let us first convince the reader that the observed nonequillibrium noise is not an artifact 
resulting simply from contacts influence or overheating of the sample by flowing current, 
as indicated above. At low temperatures the resistivity of the LCMO crystal increases by 
more than an order of magnitude. Therefore, voltage contacts shunting effects may become 
significant, in particular for contacts separated by distances comparable, or smaller than 
the sample thickness, as in our case. However, one expects shunting to manifest itself as 
a weaker than quadratic dependence of the noise power on bias, but not as a decrease of 
the noise with increasing bias, as in our experiments. The decisive test for contacts effects 
consisted in comparing the results of standard 4-point measurements with data obtained in 
an inverted contact arrangement, in which current was fed into the sample through inside 
pair of voltage contact and voltage noise measured at outside current contacts. The bias 
at which the onset of nonequillibrium noise appears in the inverted contact arrangement 
was found to be the same as the onset bias measured using standard contact arrangement. 
Moreover, in early experiments we have used much bigger, 2.6 mm separation between the 
voltage contacts. The measured value of low temperature resistivity of our sample was not 
dependent on the distance between the voltage contacts. Together, the above observations 
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enabled us to reject the possibility that low temperature noise behavior is associated with 
contact artifacts. 

The increase of the resistance by more than an order of magnitude leads to increase of the 
the power dissipated in the sample with respect to the power dissipated at higher tempera- 
tures at the same current flow. This is clearly seen in Fig. H] where the nonequillibrium noise 
at 79 K appears at bias voltages, and power levels, significantly exceeding those reached 
at higher temperatures. Temperature independence of the normalized spectral density of 
resistivity fluctuations S r causes the measured noise to follow the temperature evolution of 
the sample resistance squared. At temperatures below « 100 K the resistance decreases with 
increasing temperature, see Fig. [TJ Therefore, one may suspect that decrease of the noise 
at 79 K for currents above 1 mA and deviations from the equilibrium quadratic scaling may 
be due to simple overheating by current flow. Indeed, it is possible to explain qualitatively 
the noise behavior if one assumes that current flow above 1 mA increases the temperature 
of the sample, thus causing the resistivity decrease, what should lead to the decrease of the 
noise. This scenario requires however huge overheating of the sample by almost 100 K. This 
is unlikely in our experimental arrangement. Moreover, the changes in sample resistance 
measured at 79 K as a function of current are negligible with respect to those required to 
explain the noise behavior by overheating. Therefore, a simple overheating scenario outlined 
above cannot explain our data. 

Nevertheless, to completely reject the plausibility of the overheating scenario we have 
performed independent measurements of the sample magnetic susceptibility under current 
flow. Figure H shows the real part of magnetic susceptibility as a function of dc current, while 
the insert demonstrates the temperature dependence of susceptibility measured at different 
currents. It is evident that at T =77 K there is noticeable increase of the sample temperature 
for currents above 1 = 8 mA. This current level is almost an order of magnitude higher 
than the threshold current at which nonequilibrium noise appears. Therefore, independent 
measurements of temperature and current dependence of magnetic susceptibility enable us 
to fully reject the overheating scenario. 

To find the real origin of the noise evolution with increasing bias at 79 K, we have 
performed the measurements of the voltage dependence of the differential conductance G(V) 
and compared it with Eq. [TOJ We have determined the conductance exponent n by fitting 
the G(V) characteristics to Eq. [TD] in several voltage ranges and verified the results by 
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FIG. 7: Color online. Real part of magnetic susceptibility as a function of the current bias. Inset 
shows temperature dependence of susceptibility with 7=0, 0.1, and 10 mA current flow in the 
sample. 

differentiating the experimental curves G(V) and calculating the value of n as 

Figure |8] shows details of the dependence of 1/ f noise intensity and sample conductance 
on bias voltage at 79 K, together with the evaluated values of exponents n. We find that 
for voltages below 195 mV n = 1.81 ± 0.02. Most probably, the dominating tunneling 
mechanism in this voltage range is direct elastic tunneling across the insulating interlayer 
between the conducting banks, see the previous section. For the bias range between the two 
first maxima of the noise, i.e. between 195 and 450 mV, the dependence G(V) is almost 
linear with n = 0.98 ± 0.01. It has been argued 23 that a linear conduction background with 
n = 1 appears due to strong direct inelastic tunneling involving a broad continuum density of 
states of bosonic excitations inside the tunneling barrier. This situation has been frequently 
experimentally observed in disordered perovskite materials like high-T c superconductors and 
various CMR manganites.— 

For voltages between 450 mV to 620 mV, i.e., above the second peak of the noise, we find 
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that G(V) is composed of two shifted nonlinear pieces, each with n = 1.36 ± 0.02, giving 
the overall exponent in this range n = 1.52 ± 0.02. At voltages between 620 and 870 mV, 
i.e., above the second peak, we find n = 2.5 ± 0.03. The two last voltage ranges correspond 
to indirect inelastic tunneling via more than two conducting inclusions inside the insulating 
interlayer.— 

Changes of the exponent n reflect changes in the conduction mechanism. As we have dis- 
cussed in details previously, the electron transport in Lao^Cao.isMnOs at low temperatures 
is dominated by intrinsic tunneling mechanism.— It is therefore tempting to analyze the be- 
havior of the exponent n in terms of changing tunneling mechanisms. This approach encoun- 
ters, however, formal difficulty associated with the fact that our sample does not constitute 
a well defined single tunnel junction. The tunneling conductance in bulk Lao.82Cao.i8Mn0 3 
is associated with intrinsic tunnel barriers appearing spontaneously along the ferromagnetic 
metallic percolation paths.— The exact nature and individual properties of the barriers can- 
not be characterized as one has an experimental access only to some averaged ensemble-like 
tunnel characteristic of the sample. Therefore, we can only approximate the characteristics 
of intrinsic junctions with the known characteristics of well defined discrete junctions. 

The scenario in which the noise intensity varies with changing bias due to changes in 
conduction mechanism can be outlined as follows. At low bias voltages a direct tunneling 
mechanism is dominating, n is close to 2 for dominantly elastic charge transmissions and 
close to 1 in the inelastic tunneling regime. Crossover to nonequilibrium noise is associated 
with the change of n from 2 to 1, i.e., from elastic to inelastic transmission, which occurs 
when a — 1. At higher voltages the probability of indirect inelastic tunneling is much higher 
since the number of involved localized states in the barrier grows with increasing bias. 26 The 
conductivity of an inelastic channel increases exponentially with increasing number N of the 
localized states in the barrier. Therefore, the current is shunted by channels with higher 
N whenever available. The experimentally determined exponent n — 1.5 corresponds to 
admixture of indirect tunneling through channels with two (n = 1.3) and three (n = 2.5) 
localized states in the barrier. A similar mechanism has been invoked previously to explain 
unusual decrease of the noise with increasing bias in magnetic tunnel junctions. 15 

Figure [9] shows how the spectral exponent a at 79 K behaves with increasing voltage bias. 
The exponent reflects directly changes in the derivative d(S(f,T)f)/dlnf obtained from 
the slope of S'(/)/-vs-ln f / f characteristics and related to a by Eq. (jlj). Figures [S] and M 
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FIG. 8: Color online. Effect of the voltage bias on the differential conductance (a solid line) and 
the normalized noise spectrum at frequency 1 Hz (circles) recorded at 79 K. Vertical lines show 
voltage intervals where the exponent n in Eq. [11] are nearly constant, corresponding values of n 
are indicated in the figure. 

evidence apparent correlations between the behavior of the noise, evolution of the exponent 
a, and corresponding n values. It seems that significant changes in the difference 5a = a — 1, 
or equivalently the changes of the sign of the derivative d(S(f,T)f)/d\nf, appear in close 
vicinity of voltages for which 5a = 0, i.e. for voltages at which the distribution function 
-P(A) is a constant independent on the magnitude of an activation energy. 

Changes of the conduction mechanism, reflected by changes in n, are accompanied by 
clear noise peaks. The peaks can be seen as manifestation of excess partition noise at specific 
voltages at which alternative inelastic channels with different iV are available for electron 
transfer. We conclude that Figs. [8] and [9] evidence a clear correlation between bias induced 
changes in the noise and modifications of / — V characteristics reflecting changes in the 
intrinsic tunneling mechanism. 

The explanation of noise decrease with increasing bias can be straightforwardly based 
on a scenario in which bias causes changes in the symmetry of two level ffuctuators in the 
ensemble responsible for generating the observed noise. Elementary ffuctuators in intrinsic 
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FIG. 9: Color online. Effect of the voltage bias on the parameter a at frequency 1 Hz (circles) 
recorded at 79 K. Vertical lines show voltage intervals where the exponent n in Eq. [11] are nearly 
constant, corresponding values of n are indicated in the figure. 

tunnel junctions can be associated, for example, with charge traps within tunnel barriers and 
two distinct states of a fluctuator with empty and loaded state of the trap, or alternatively, 
with charged oxygen vacancies jumping between two sites. Electrostatic field of the charge 
trap modulates the height of the tunnel energy barrier and, consequently, the resistivity of 
the junction, according to the trap occupancy state. 

An alternative mechanism consists in bias induced changes to the regions where relax- 
ation of tunneling charges takes place. At low bias, when elastic tunneling dominates, the 
transmitted electrons relax energy in conducting regions through interactions with bosonic 
degrees of freedom, whereas for inelastic tunneling the relaxation occurs inside the insulating 
barriers. If in the latter case the activation energy of defects within the barrier is higher 
than that of the conducting regions, then the noise will decrease with increasing intensity 
of the inelastic tunneling events. The observed bias induced changes in -P(A) distribution 
with changing bias do allow, at this stage, to choose one from the above discussed scenarios. 
One has however to note that the energy threshold corresponding to the onset on nonequlib- 
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rium noise at threshold bias of 200 mV is much lower than the activation energy of oxygen 
vacancies in perovskites. At the same time we observe that this is the energy level of the 
Jahn- Teller distortions in manganites. 32 

V. CONCLUSIONS 

In conclusion, we have investigated conductivity fluctuations in current biased 
La .82Ca .i8MnO 3 single crystals in a wide range of bias current and temperature, a very 
good model object for studying the noise processes in inhomogeneous conducting media since 
different magnetic and electric transport properties within a temperature interval from the 
liquid-nitrogen temperature up to ambient one. The observed noise retains 1/f spectra in 
entire experimentally explored range of bias and temperature despite pronounced changes 
in the dissipation mechanism and magnetic state that occur in the same range. The voltage 
noise measured at constant current bias J was found to have characteristics of a quasi- 
equilibrium modulation noise with intensity scaling as I 2 , indicating that the conductivity 
fluctuations are only revealed by the current flow and not stimulated or modified by it. Only 
at liquid nitrogen temperatures the observed 1/f noise behaves like nonequilibrium one and 
its intensity not only does not scale with the square of the bias but even decreases locally 
with increasing bias. 

We have verified that in entire experimentally investigated range the noise fulfils the 
Dutta-Horn model frequency-temperature reciprocity condition. We have concluded there- 
fore that noise in our system arises from a superposition of thermally activated processes 
with broad distribution of activation energies. However, basing this conclusion just on the 
consistency between temperature and frequency dependence of the noise has to be taken 
with the grain of salt since the noise model based on serious kinetics of a random walk in 
random potential (RWRP), a model that in principle could apply well to our 0.18 LCMO 
sample, produces frequency-temperature reciprocity relation identical to Eq. (JSJ). Only in- 
depth analysis of non-Gaussian components of the noise, provided they exist, enables one to 
distinguish between the reality of two models.- We did not observe any non-Gaussian noise 
components in the experiment and are unable to distinguish between the two models. 

Nevertheless, we have tentatively accepted the conclusion that the noise indeed arises 
from the DH mechanism and inferred the distributions of activation energies in the elemen- 
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tary fluctuations ensembles by analyzing the power spectral density of the noise within the 
framework of the DH approach. We have found that the appearance of a nonequilibrium 
noise is directly linked with the onset of bias dependence of the distribution of activation 
energies P{E). We have found clear correlations between the changes of the energy distri- 
bution and the power exponent of the I — V curves. Changes in the power exponent, i.e., 
the hallmarks of the changes in the dissipation mechanism, are correlated with changes in 
the noise behavior. We have discussed two possible scenarios: (i) a direct influence of the 
electric field on fluctuators which in this case should be charged and (ii) an effect of the 
spatial change of the region where the electron energy dissipates during its traveling across 
an inhomogeneous sample. Both models are appropriate in our case and we are performing 
additional studies in order to reveal the most probable one. 

The results obtained in the work open a perspective for applying the noise tool to inves- 
tigate the nature and physics involved in metastable resistivity states which can be enforced 
in low-doped manganites by applying strong pulsed electric field/current at low tempera- 
ture.— 1 ^ Our first experiments have confirmed this assumption and their analysis will be 
published elsewhere. 
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